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Abstract. Many of the high mass X-ray binaries (HMXRBs) discovered in recent years in our Galaxy are characterized by 
a high absorption, most likely intrinsic to the system, that can impede their detection at the softest X-ray energies. Exploiting 
the good coverage obtained with sensitive XMM-Newton observations, we have undertaken a search for highly absorbed X-ray 
sources in the Small Magellanic Cloud (SMC), which is known to contain a large number of HMXRBs. After a systematic 
analysis of 62 XMM-Newton SMC observations, we obtained a sample of 30 sources with evidence of an equivalent hydrogen 
column density larger than 3x10^"^ cm~^. Five of these sources are clearly identified as HMXRBs, four being previously 
known (including three X-ray pulsars) and one, XMMU J005605. 8-720012, being reported here for the first time. For the 
latter, we present optical spectroscopy confirming the association with a Be star in the SMC. The other sources in our sample 
have optical counterparts fainter than magnitude ~16 in the V band, and many have possible NIR counterparts consistent 
with highly reddened early-type stars in the SMC. While their number is broadly consistent with the expected population of 
background highly absorbed active galactic nuclei, a few of them could be HMXRBs in which an early- type companion is 
severely reddened by local material. 
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5 1. Introduction 



Many observations of the Small Magellanic Cloud (SMC) in 
the X-ray energy band have led to the discovery of a large num- 
ber of High Mass X-ray Binaries (HMXRBs). The number of 



known HMXRBs in the SMC (about one hundred, iLiu et al. 



(l2005h ) is much larger than expected by scaling the number 
of these sources seen in the Milky Way according to the mass 
ratio of the two galaxies [Mmw^^smc ~50). This has been 
interprete d as evidence of a recent episode of star formation i n 
the SMC (iMai id et al.' (200^), IShtvkovskiv & Gilfanovl (|2007'), 



Antoniou et al. (.2010.) ). It is remarkable that the SMC contains 



only one supergiant system, SMC X-1. All the other HMXRBs 
in the SMC consist of neutron stars accreting from Be type 
companions. Be X-ray binaries constitute also the largest class 
of HMXRBs in our Galaxy, but the relative number of super- 
giant systems is much h igher thaii in the SMC, since 30 of the 
1 14 Galactic HMXRBs (JLiu et al.ll2006l) are confirmed or sus- 
pected supergiant systems. 



Send offprint requests to: Giovanni Novara, novara@iasf- 
milano.inaf.it 



Many new Galactic HMXRBs have been discovered in the 
past few years with the INTEGRAL satellite, thanks to high 
sensitivity in the hard X-ray range, coupled with an exten- 
sive monitoring of the Galactic plane. These observations led 
to the recognition of the new class of Supergiant Fast X-ray 
Transients (SFXT) and to the discovery of persistent supergiant 
systems characterized by a high absorption (equivalent column 
density above a few 10^'^ cm^^), which escaped an earlier dis- 
covery because the high absorption severely suppresses their 
flux below lOkeV (see, e.g.. Sidoli (2010) . for a recent review). 



To search for highly absorbed HMXRBs in the SMC we 
carried out a dedicated analysis of the XMM-Newton data 
collected during the large p r ograrn for the SMC survey per- 
formed by iHaberl & PietschI (120081) and other observations of 
SMC targets. In Sect. |2l we describe the X-ray observations 
and data processing and in Sect. [3] we describe the source de- 
tection and selection criteria that led us to identify a sample 
of 30 highly absorbed sources. Five of them are confirmed 
HMXRBs, including four that were already known as X-ray 
binaries (Sect. |4|i. In Sect. |5] we present optical spectroscopy 
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of the new source, demonstrating that it is a Be system in the ^ 

SMC. We finally discuss our results in Sect. |6] 

2. X-ray observations and data processing 

We analyzed the 62 XMM-Newton pointings of the SMC 
performed from 2001 May 31 to 2010 March 16 listed in 
Table [T] The three EPIC focal plane cameras (ITurner et alJ 
200 It IStriider et alj l200lh were operated in standard Full 
Frame mode. 

The data of each camera were processed independently 
using the standard XMM-Newton Science Analysis Software 
{SAS, v.10.0.0). We first used the tasks emproc and epproc 
to obtain the MOS and pn event files of each pointing. We then 
filtered out time intervals affected by high instrumental back- 
ground induced by flares of soft protons (with energies below 
a few hundred keV). To avoid contributions from genuine X- 
ray source variability, this was done by examining light curves 
binned at 100 s for events with energies above 10 keV and 
PATTERN <4 or for, respectively, the MOS peripheral CCDs 
and the whole pn camera; the good time intervals (GTI) were 
selected by adopting different count-rate thresholds for the 
different pointings and instruments. By selecting only events 
within GTI, we finally obtained three "clean" event lists for 
each observation (i.e. two for the MOS and one for the pn cam- 
eras), whose exposure times are reported in Table [T] The total 
net exposure times are respectively 1572 and 1598 ks for the 
MOS cameras, and 1361 ks for the pn one. 



3. Source detection and selection criteria 

We optimized our source detection strategy to select highly ab- 
sorbed X-ray sources. To this aim, we simulated the expected 
EPIC count distributions for absorbed power-law spectra, with 
photon-index F in the range 1-2. We used two absorption com- 
ponents: the first one, with column d ensity fixed at 6 x 10^" 
cm~^ and elemental abundances from I Wilms et al.l (120001) . ac- 
counts for the foreground absorption in our Galaxy; for the sec- 
ond one, which accounts for the absorption in the SMC and is 
intrinsic to the source, we considered Nh values in the range 
10^^-10^^ cm~^, and metal abund ances of 0.2, as is typical of 
the SMC dRussell & Dopitalll992 l). We defined the Hardness 
Ratio HR ={H - S)/{H + S)), based on the count rates in the 
soft (S, 1-3 keV) and hard (H, 3-10 keV) energy ranges. These 
simulations showed that, for F = 2 and Nh > 3 x 10^'^ cm^^, 
the fraction of the detected source counts with energies above 
3 keV is at least 70 %, and that this fraction is even higher 
for lower values of F. In Fig. [T] we report the estimated HR 
values as a function of Nh, for both the pn and MOS cameras 
and different values of F. On the basis of these simulations, we 
adopted a HR threshold value of 0.7, in order to select sources 
withNn > 3 X lO^^ cm^^. 

For each observation in Table [T] we used the cleaned event 
file to produce merged images from the three cameras in the 
1-3 and 3-10 keV e nergy range s , acco rding to the procedure 
described in detail in iBaldi et all (120021) . For the 3-10 keV pn 
images, we excluded the energy range 7.8-8.2 keV, in order 
to reduce the instrumental background caused by the Cu line. 




Fig. 1. Expected hardness ratio (of the count rates in the energy 
ranges 1-3 to 3-10 keV) as a function of the column density. 
Solid and dashed lines refer to the pn and the MOS cameras, re- 
spectively. Red, green, and blue lines refer to power-law spectra 
with photon indices of 2, 1.5, and 1, respectively. The horizon- 
tal line at HR = 0.7 represents the threshold value used to select 
highly absorbed sources. 

For each image, we also produced the corresponding merged 
exposure map, which accounts for mirror vignetting and ef- 
fective field of view. The total count-rate-to-flux conversion 
factors were obtained as a mean of the M0Sl/M0S2/pn fac- 
tors, weighted on the individual exposures of the three cam- 
eras. The source detection was based o n a maximum likeli - 
hood technique, as described in detail in iNovara et al.l (l2006h . 
and used backgrou nd maps produced with the correction algo- 
rithm described in iBaldi et alj (12002!) . For each energy band, 
we only selected sources with a final detection likelihood -In P 
> 8.5 (where P is the probability of a false detection due to a 
Poissonian random fluctuation of the background). This likeli- 
hood threshold corresponds to a '^3 a detection. After manu- 
ally removing a few spurious detections caused, for instance, 
by events falling close to the CCD edges, we obtained a fi- 
nal master list containing ~1500 sources. For each source, the 
master list provides various parameters including the detec- 
tor and sky coordinates, the effective exposure time, the total 
counts, count-rate and errors in the different energy ranges, 
and the detection likelihood. 
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Table 1. Log of the XMM-Newton observations of 62 SMC fields with the corresponding net good time interval (GTI) for the 
two MOS and the pn cameras. 





Observation ID 


XMM-Newton 


Date 




GXI (ks) 






revolution 


(UT) 


MOSl 


M0S2 


pn 


0011450101 


0270 


2001-05-31X02:20:10 


42.8 


43.0 


40.8 


0011450201 


0355 


2001-11-16X03:23:32 


37.2 


38.5 


40.1 


0018540101 


0357 


2001-11-20X23:42:37 


26.7 


27.0 


23.4 


0084200101 


0422 


2002-03-30X13:48:28 


8.8 


8.2 


7.1 


0084200801 


0340 


2001-10-17X10:07:40 


20.7 


20.7 


17.1 


0110000101 


0156 


2000-10-15X15:18:28 


25.2 


25.7 


14.7 


0110000301 


0157 


2000-10-17X21:45:54 


12.6 


12.4 


7.4 


0112780201 


0143 


2000-09-19X02:05:37 


3.6 


2.3 


0.0 


0112780601 


0254 


2001-04-29X21:07:59 


0.9 


1.5 


2.9 


0135721701 


0721 


2003-11-16X06:12:02 


24.4 


24.8 


29.8 


0142660801 


0721 


2003-11-17X03:55:54 


7.2 


7.5 


6.8 


0157960201 


0737 


2003-12-18X14:32:45 


17.9 


18.7 


17.2 


0164560401 


0803 


2004-04-28X22:08:04 


0.2 


0.1 


0.1 


0301170101 


1151 


2006-03-22X21:39:54 


18.8 


19.7 


16.4 


0301170201 


1151 


2006-03-23X04:48:17 


22.5 


22.7 


14.4 


0301170301 


1158 


2006-04-06X04:32:35 


17.9 


18.0 


14.4 


0301170601 


1153 


2006-03-27X12:21:01 


14.4 


14.5 


11.7 


0311590601 


1146 


2006-03-13X15:17:13 


11.2 


11.0 


9.6 


0402000101 


1248 


2006-10-03X00:09:09 


21.7 


21.7 


19.0 


0403970301 


1329 


2007-03-12X20:02:20 


28.0 


28.7 


19.2 


0404680201 


1263 


2006-11-01X00:56:29 


32.4 


32.4 


30.8 


0404680301 


1344 


2007-04-11X19:38:25 


20.1 


19.9 


15.4 


0404680501 


1344 


2007-04-12X03:07:23 


23.3 


23.7 


21.6 


0500980101 


1380 


2007-06-23X05:51:39 


25.2 


25.0 


22.3 


0500980201 


1372 


2007-06-06X08:52:16 


27.6 


28.3 


14.1 


0501470101 


1371 


2007-06-04X08:59:50 


13.4 


16.5 


9.6 


0503000201 


1444 


2007-10-28X05:49:58 


21.5 


21.5 


20.0 


0503000301 


1514 


2008-03-16X15:25:16 


17.9 


19.3 


2.0 


0601210101 


1727 


2009-05-14X09:14:24 


20.7 


21.2 


19.8 


0601210201 


1794 


2009-09-25X00:15:42 


37.5 


37.5 


35.7 


0601210301 


1729 


2009-05-18X10:29:15 


30.9 


31.4 


27.5 





Table 1. Continued 
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Observation ID 


XMM-Newton 


Date 




GXI (ks) 






revolution 


(UT) 


MOSl 


M0S2 


pn 


0601210401 


1794 


2009-09-25X11:22:24 


37.6 


37.6 


36.0 


0601210501 


1794 


2009-09-25X22:30:43 


49.8 


50.4 


39.9 


0601210601 


1795 


2009-09-27X00:10:53 


36.7 


37.1 


32.7 


0601210701 


1795 


2009-09-27X11:35:54 


38.7 


38.6 


37.1 


0601210801 


1801 


2009-10-09X18:33:30 


24.7 


24.7 


23.1 


0601210901 


1795 


2009-09-27X23:02:34 


34.9 


35.3 


33.9 


0601211001 


1817 


2009-11-09X21:16:08 


34.0 


35.0 


23.7 


0601211101 


1806 


2009-10-18X22:47:09 


31.6 


31.6 


30.1 


0601211201 


1807 


2009-10-20X22:51:09 


33.1 


33.2 


29.0 


0601211301 


1798 


2009-10-03X05:08:47 


32.4 


32.5 


30.9 


0601211401 


1814 


2009-11-04X21:38:31 


35.3 


35.6 


31.1 


0601211501 


1803 


2009-10-13X00:02:01 


37.6 


37.4 


34.0 


0601211601 


1802 


2009-10-11X22:43:47 


31.6 


32.3 


28.8 


0601211701 


1804 


2009-10-16X01:05:21 


25.1 


26.6 


18.4 


0601211801 


1819 


2009-11-13X20:59:56 


26.2 


30.2 


22.3 


0601211901 


1827 


2009-11-30X14:46:19 


31.7 


31.7 


30.1 


0601212001 


1826 


2009-11-27X22:39:06 


28.3 


28.3 


26.3 


0601212101 


1820 


2009-11-16X05:48:08 


34.1 


34.1 


32.5 


0601212201 


1822 


2009-11-16X06:10:27 


27.1 


29.4 


24.3 


0601212301 


1786 


2009-09-09X09:13:53 


33.5 


33.5 


31.9 


0601212401 


1750 


2009-06-29X14:46:19 


29.5 


30.9 


23.9 


0601212501 


1786 


2009-09-09X19:13:54 


33.5 


33.5 


31.9 


0601212601 


1750 


2009-06-29X06:04:39 


26.0 


27.1 


17.5 


0601212701 


1840 


2009-12-26X07:25:22 


36.7 


36.7 


34.5 


0601212801 


1831 


2009-12-07X23:35:54 


21.3 


21.3 


16.7 


0601212901 


1788 


2009-09-13X13:29:26 


36.1 


36.1 


34.5 


0601213001 


1788 


2009-09-13X01:11:03 


41.6 


41.7 


38.9 


0601213201 


1878 


2010-03-12X00:56:15 


12.5 


13.2 


9.4 


0601213301 


1878 


2010-03-12X05:26:15 


9.8 


9.9 


8.1 


0601213401 


1880 


2010-03-16X10:05:12 


17.6 


17.5 


11.9 


0656780201 


1886 


2010-03-27X12:20:41 


11.0 


12.3 


7.0 





Table 2. Main characteristics of the highly absorbed X-ray sources. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


SRC 


NAME 


RA 


DEC 


ERROR 


Likelihood 


RATE 


fx 


HR 


Remarks 






h:m:s 


o:/:// 


arcsec 




(xlO^'^ctss^^) 


erg cm~^ s~^ 







01 


XMMUJ004102.6-732530 


00:41:02.69 


-73:25:30.9 


2.1 


27.43 


02 


XMMUJ004210.9-73301 1 


00:42:10.90 


-73:30:11.8 


2.0 


16.80 


03 


XMMUJ004226.3-730417 


00:42:26.38 


-73:04:17.4 


1.5 


37.07 


04 


XMMUJ004244.7-732356 


00:42:44.77 


-73:23:56.4 


2.1 


15.34 


05 


XMMUJ004536.2-724131 


00:45:36.23 


-72:41:31.3 


1.7 


10.74 


06 


XMMUJ0048 14.0-73 1006 


00:48:14.07 


-73:10:06.0 


1.5 


58.28 


07 


XMMUJ004818.7-732102 


00:48:18.74 


-73:21:02.4 


1.9 


>86 


08 


XMMUJ004853.5-732455 


00:48:53.51 


-73:24:55.7 


2.0 


>86 


09 


XMMUJ00491 1.5-731717 


00:49:11.57 


-73:17:17.8 


1.8 


12.79 


10 


XMMUJ005020.7-720907 


00:50:20.74 


-72:09:07.0 


0.6 


>86 


11 


XMMUJ005306.6-722400 


00:53:06.66 


-72:24:00.2 


1.5 


12.96 


12 


XMMUJ005322.4-7 15927 


00:53:22.44 


-71:59:27.5 


1.0 


10.35 


13 


XMMU J005 605 . 8-7200 1 2 


00:56:05.85 


-72:00:12.2 


1.6 


26.57 


14 


XMMUJ005722.4-7 13114 


00:57:22.44 


-71:31:14.8 


0.9 


12.66 


15 


XMMUJ005724.3-715917 


00:57:24.31 


-71:59:17.8 


1.5 


16.81 


16 


XMMUJ005732.5-712926 


00:57:32.51 


-71:29:26.3 


1.8 


8.93 


17 


XMMUJ005735.8-721935 


00:57:35.82 


-72:19:35.1 


2.0 


30.99 


18 


XMMUJ0058 12.9-723049 


00 


58:12.92 


-72:30:49.5 


1.4 


>86 


19 


XMMUJOlOl 15.2-721640 


01 


01:15.28 


-72:16:40.6 


2.1 


13.14 


20 


XMMUJ010233.9-723443 


01 


02:33.99 


-72:34:43.2 


1.8 


60.86 


21 


XMMUJ010248.6-730822 


01 


02:48.69 


-73:08:22.6 


1.4 


39.22 


21 


XMMUJ010432.5-722543 


01 


04:32.56 


-72:25:43.6 


2.1 


12.73 


23 


XMMUJ010802.9-722627 


01 


08:02.93 


-72:26:27.5 


2.1 


17.50 


24 


XMMUJ0108 1 1 .9-721005 


01 


08:11.91 


-72:10:05.2 


1.2 


24.42 


25 


XMMUJO 1083 1 .4-730630 


01 


08:31.40 


-73:06:30.9 


1.6 


39.61 


26 


XMMUJ010842.6-723839 


01 


08:42.68 


-72:38:39.6 


1.8 


25.94 


27 


XMMUJ010908.0-720642 


01 


09:08.09 


-72:06:42.1 


1.3 


13.52 


28 


XMMUJOl 1339.7-724733 


01 


13:39.74 


-72:47:33.3 


2.0 


11.87 


29 


XMMUJOl 1434.8-730730 


01 


14:34.83 


-73:07:30.4 


2.2 


10.66 


30 


XMMUJOl 1807.8-73 1713 


01 


18:07.84 


-73:17:13.9 


2.2 


16.17 



13.6 ±2.9 
12.4 ± 3.3 

11.7 ±2.0 
5.9 ± 1.6 

18.4 ±5.8 

55.4 ±8.2 

103 ±9.1 

84.3 ± 8.2 

15.1 ±4.8 

68.6 ± 5.3 
6.0 ± 1.8 

11.2 ±3.6 
8.4 ± 1.7 

15.7 ± 4.3 

17.8 ±4.6 

11.3 ±3.7 

21.4 ±4.5 

1330 ±30 

19.7 ±5.5 

29.7 ± 4.9 

28.2 ± 6.2 

8.9 ± 2.5 

7.9 ± 2.2 

7.0 ± 1.6 

12.1 ±2.1 

39.5 ± 8.0 
6.3 ±1.9 

15.0 ±4.8 
14.7 ±4.7 

24.6 ± 6.3 



8.94E-14 
8.21E-14 
7.76E-14 
3.89E-14 
1.18E-13 
3.49E-13 
6.81E-13 
5.71E-13 
9.86E-14 
4.51E-13 
3.93E-14 
7.41E-14 
5.51E-14 
1.03E-13 
1.20E-13 
7.48E-14 
1.45E-13 

8.78E-12 
1.30E-13 
1.90E-13 
1.83E-13 
5.81E-14 
5.21E-14 
4.61E-14 
7.98E-14 
2.56E-13 
4.11E-14 
l.OlE-13 
9.70E-14 
1.60E-13 



0.86 ± 0.23 
1.00 ±0.29 
0.73 ± 0.20 
0.80 ± 0.34 
0.73 ± 0.29 
0.73 ±0.12 
0.81 ±0.08 
0.72 ± 0.07 
0.91 ± 0.34 
0.94 ± 0.05 
0.72 ± 0.33 
0.87 ± 0.36 
0.70 ± 0.23 
0.72 ± 0.33 
0.84 ± 0.29 
0.86 ± 0.46 
1.00 ±0.18 

0.77 ± 0.01 
0.82 ± 0.35 
0.91 ±0.17 
0.77 ± 0.23 
0.85 ± 0.36 
0.71 ±0.30 
0.81 ±0.27 
0.93 ±0.19 
0.79 ±0.17 
0.85 ± 0.35 
0.86 ±0.31 
0.70 ± 0.29 
0.85 ± 0.23 



XMMUJ004814.1-731003, SXP25.5 (1) 
[SG2005] SMC 34 (2) 

CXOU J004910.7-731717 (3) 



New Be binary 



XMMU J005735.7-721932 (4) 
CXOU J005736.2-72 1934, SXP565 (5) 
RX J0058.2-7231, SXP293 (6) 



P 



1: 



X 



a' 



in 
O 



Key to Table: Col.(l) = source ID; Col. (2) = source catalogue name; Col.(3) = right ascension (J2000); Col. (4) = declination (J2000); Col. (5) = position error (1 a c.l.); Col. (6) = detection 
likelihood; Col. (7) = 3-10 keV count rate (sum of the FN and MOS counts divided by the sum of FN and MOS exposures); Col. (8) = X-ray flux in the energy band 3-10 keV, assuming a 
power-law emission model with a hydrogen column density Nh > 3 x 10^'^ cm~^ and a photon-index F = 1.5; Col.(9) = hardness ratio between energy bands 1-3 keV and 3-10 keV; Col.(lO) 
= lik ely identific a tion wi th al ready known sources. Reference s are given in parenthesis: 

(1) iHaberl et alJ (l2008bl) (2) IShtvkovskiv & Gilfanovl ( |2005|) (3) JLavcock et ai] ( 1201 Oh (4) ISasaki et alj ( |2003|) , iHaberl et al.l l l2008ah (5) iMacomb et al] ( l2003h (6) iHaberl et al.l ( l2008ah (7) 
iMcGowan eTarf l l2008l) 
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4. Results 




Hardness Ratio 

Fig. 2. Hardness-ratio distribution of the detected X-ray 
sources 

In Fig. 121 we report the HR distribution of all the detected 
sources. This figure illustrates that most sources have HR < 0, 
with a peak around -0.5, hence they are characterized by non- 
absorbed soft spectra. On the other hand, only very few sources 
have high HR values. In our final source list, there are 30 ob- 
jects with HR > 0.7, which we consider as highly absorbed 
sources with Nh > 3 x 10^^ cm~^; given the HR distribution, 
we expect that, even if we take into account the HR errors, only 
a very few additional sources would have HR > 0.7. The sky 
coordinates of these 30 objects, corrected for the satellite point- 
ing uncertainty, are listed in Tabled The astrometric correction 
of the X-ray images was done by cross-matching the brightest 
sources of each field with the Ma gellanic Clouds Photo metric 
Survey optical catalogue (MCPS. IZaritskv et all (l2002h ), and 
selecting the X-ray sources with a single, clearly evident opti- 
cal counterpart within a 3" error radius. We then used the IRAF 
task geomap to compute the linear transformation between the 
X-ray and optical coordinates of these reference sources. The 
calculated X-ray-to-optical coordinate transformations have 
rms fit residuals between 0'.'43 and 1'.'85 at 1 cr c.l., with an av- 
erage value of 1'.' 1 ; they correspond to the (radial) systematic 
astrometric errors in the corrected observations. These trans- 
formations were then used by the IRAF task geoxytran to 
obtain the corrected X-ray coordinates given in Table |2l The 
total uncertainty on the X-ray position of each source is given 
by the sum in quadrature of the calculated systematic astromet- 
ric error in the observation and the measured statistical error of 
the source itself. 

We also list in Table |2] the hardness ratios, count rates, and 
fluxes of the sources. The fluxes refer to the 3-10 keV energy 
range, and have been computed assuming a power-law with 
photon index F = 1.5 and Nh = 3x10^"' cm^^. The brightest 
source has a absorbed flux of 9x10^^^ erg cm^^ s^^, while 
the faintest sources in our sample have fluxes of ^4x10^^'* 
erg cm~^ s~^ correspon ding to ~2 x 10'^'^ erg s~^, assuming a 



SMC distance of 60 kpc JHilditch et al.ll2005l) 



Some of our selected sources were already reported in other 
studies of the SMC, in these cases we give in the table their 
original name and the relevant references. Four sources, includ- 
ing three pulsars, were already known as HMXRBs: SXP25.5, 
SMC34, SXP565, and SXP293. We searched for possible coun- 
terparts of our highly absorbed sources using the MCPS cata- 
logue, which is based on observations performed in the U, B, 
y, and / filters betw een November 1996 and December 1999 
JZaritskv et al.ll2002h . and the InfraRed Survey Facility (IRSF) 
MageUanic Cloud catalogue (.Kato et aL,2007,) for the NIR (7, 
//, and /i: filters). 

For each source, we report in Table [3] the optical and NIR 
properties of the brightest object in the K band present within 
a distance equal to two times the 1 a position error (with only 
five exceptions, marked in the table, the brightest object in K is 
also the brightest in V). For five sources, no catalogued optical 
or NIR counterparts were found to limits as faint as of V~23 
mag and K~17 mag. 

A plot of the NIR versus optical colors shows that the four 
already known HMXRBs of our sample are located in the re- 
gion corresponding to J-K<0.5 mag and B-V<0.5 mag (see 
Fig. O. As discussed in Sect. 5, the other source falling in 
this region (source # 13) is a new Be binary system in the 
SMC. Three of these sources are well above the po sition of 
the ma in-sequence stars (whose colors are taken from 'johnsonl 
(11966 )): if they are affected by an IR excess (which could be 
due to the decretion disk of the Be star), their optical/NIR col- 
ors are consistent with early-type stars with moderate optical 
reddenings (Ay ^ 2-3 mag). On the other hand, HMXRBs 
in which the optical companion is also affected by a high lo- 
cal absorption would show much redder colors, as seen in sev- 
eral of the new HMXRBs discovered by INTEGRAL in our 
Galaxy. Therefore, we cannot exclude that some of the other 
sources plotted in Fig. [3] or sources for which no optical and 
NIR counterparts have been detected, are also highly absorbed 
HMXRBs. The five sources at the top of the plot are far from 
the main-sequence stars, therefore it is difficult to suggest a 
reliable classification of them; this is particularly true for the 
three sources shown as open circles, for which a different coun- 
terpart could also be suggested. These five sources could be 
early-type stars only if they are characterized by both a large 
IR excess and a hi gh reddening. On the other hand, since for 
AGNs J- K > .5 (Kouz uma & Yamaokall20l d) and -0.4 < B- 
V < 0.7 dHatzi minaoglo u et al.ll2002l) . these sources could be 
more probably classified as reddened AGNs. 

To search for possible periodicities in any optical coun- 
terparts close to the XMM-Newton sources, data from OGLE 
III were acquired for all objects within 4" of the 30 XMM- 
Newton position. We applied a Lomb-Scargle analysis to a to- 
tal of 129 lightcurves, searching for periods in the range 2-200 
days. Periods shorter than 2 days were avoided because they 
are getting close to the average daily sampling of the OGLE 
data, and periods longer than 200 days approach the annual 
sampling patterns. In each case, the data were first de-trended 
using a polynomial of order three before being searched. In 
some cases, there were many optical counterparts; conversely, 
other sources had no optical objects within the search zone. 
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Table 3. Candidate optical-IR counterparts of the highly absorbed X-ray sources. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


ID 


NAME 


V 


fx/fv 


K 


B-V 


J-K 






mag 


logio 


mag 


mag 


mag 



01 


IRSF 00410291-7325332 


20.79'"^ 


0.76 


16.68 


1.21 


1.61 


02 


IRSF 00421082-7330127 


19.31<") 


0.13 


17.00 


1.05 


0.63 


03 


IRSF 00422638-7304184 


>23 


>1.58 


15.65 


- 


2.02 


04 


IRSF 00424550-7323553 


17.05 


-1.09 


13.90 


1.39 


0.87 


05 


IRSF 00453646-7241343 


19.43 


0.34 


15.72 


0.73 


1.84 


06 


IRSF 00481410-7310040 


15.30 


-0.85 


15.64 


0.26 


0.04 


07 


IRSF 00481871-7321000 


16.18 


-0.21 


15.29 


0.25 


0.36 


08 


IRSF 00485330-7324574 


>23("' 


>2.44 


15.55 


- 


1.30 


09 


IRSF 00491089-7317172 


18.80 


0.01 


15.95 


1.00 


0.38 


10 


IRSF 00502062-7209073 


>23 


2.34 


16.86 


- 


1.65 


11 


? 


>23 


>1.28 


>17 


- 


- 


12 


MCPS 2493874 


20.74 


0.65 


>17 


0.47 


- 


13 


IRSF 00560575-72001 18 


16.71 


-1.08 


17.18 


-0.11 


-0.11 


14 


7 


>23 


>1.70 


>17 


- 


- 


15 


MCPS 3076421 


21.42 


1.14 


>17 


0.35 


- 


16 


7 


>23 


>1.56 


>17 


- 


- 


17 


IRSF 00573600-7219342 


15.99 


-0.95 


15.35 


0.01 


0.42 


18 


IRSF 00581258-7230488 


14.87 


0.38 


14.14 


0.10 


0.39 


19 


IRSF 0101 1529-7216374 


19.23'"' 


0.30 


14.69 


1.28 


1.59 


20 


IRSF 01023342-7234421 


19.64'"' 


0.62 


16.09 


1.14 


1.68 


21 


IRSF 01024839-7308213 


20.99 


1.15 


17.08 


1.04 


1.90 


22 


IRSF 01043305-7225407 


16.23 


-1.25 


13.12 


1.38 


0.81 


23 


MCPS 441 1974 


20.01 


0.21 


>17 


0.04 


- 


24 


IRSF 01081 153-7210048 


>23 


>1.35 


17.26 


- 


1.61 


25 


IRSF 01083147-7306325 


18.00 


-0.41 


15.84 


0.98 


0.50 


26 


IRSF 01084310-7238375 


19.53 


0.71 


17.84 


0.84 


0.21 


27 


IRSF 01090778-7206438 


21.13 


0.55 


>17 


0.57 


- 


28 


MCPS 4879163 


22.36 


1.44 


>17 


0.56 


- 


29 


7 


>23 


>1.67 


>17 


- 


- 


30 


7 


>23 


>1.89 


>17 


- 


- 



Key to Table: Col.(l) = source ID; Col. (2) = name of the brightest infrared (IRFS sources, (Kato et al. 2007)) or optical (MCPS sources, 

jZaritskv et al.ll2002h ) candidate counterpart within the X-ray error circle (a '?' is reported when no counterpart is found); Col. (3) = V 

magnitude; Col.(4) = X-ray-to-optical flux ratio (in logio units); Col. (5) = K magnitude; Col.(6) = B-V color index; Col. (7) = J-K color index 

Notes '"' for this object there is a possible counterpart with a slightly brighter V magnitude. 



Significant periodicities were found for four sources (three 5. XMMU J005605.8-72001 2: A new Be HMXRB 

of which were already known) and tentative periods for an in the SMC 

additional three objects (see Table [Hi. That three previously 

known periods were "re-discovered" in this search reassures Source #13 was detected with a total of 68 counts in the hard 

us that the techniques used are effective in identifying possible energy range 3-10 keV and a HR value of 0.71, just above our 

periods. The one new significant detection (source number 19) threshold. It has a relatively bright (V = 16.71 mag) optical 

shows a modulation that is extremely sinusoidal in nature when counterpart, with color B-V = -0. 12 mag, and with a quality 

the data are folded at the period of 199.6 days. This modula- flag in the MCPS indicating a successful fit with a stellar atmo- 

tion shape is not normal for a long-period binary modulation, sphere model. It has an unabsorbed luminosity, in the energy 

which tends to reveal short outbursts coincident wi th the perias- range 3-lOkeV, of 4.4x10^'' erg s^^. 



tronpassage of the neutron star (see, for example,| Coe & Edge| Optical spectral observations of source # 1 3 were obtained 



(12004D). Such smooth sinusoidal modulations tend to be associ- ;„ September 2010, using the 1.9-meter telescope and the 



ated with pulsations in stars. A likely possibility is that it could Cassegrain spectrograph at the South African Astronomical 
be an Non-Radial Pulsations (NRP) in a B-type star with a pe- observatory (SAAO) in Sutherland. We used grating number 7 
riod very close to the one day sampling, with the result that (ggg Unes per mm) to obtain spectra between 3700 A and 7700 
the signal at 199.6 days is the beat frequency of the true pulse ^ ^^ ^ resolution of 5 A. For these, the slit size was 1'.'5 x 1 .'5. 
period and the sampling period. Exposure times were limited to 900 s with a positional accuracy 

of < 1". Data reduction included bias subtraction and flat-field 
correction using the IRAF software package. Extraction (task 
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Table 4. List of OGLE III objects within 4" of an XMM-Newton position that show evidence of a periodic signal 



XMM ID OGLE III ID 



V 



I LS power Period 



Comments 



2 


22605 


21.3 


20.3 


12 


2.16 days 


6 


50768 


15.8 


15.7 


14 


22.5 days 


8 


20669 


7 


21.0 


13 


4.08 days 


17 


27992 


16.1 


15.8 


26 


157.2 days 


IB 


19191 


14.9 


14.6 


200 


59.7 days 


20 


43690 


20.8 


21.0 


15 


5.96 days 


26 


22676 


19.6 


18.6 


23 


199.6 days 



Tentative 

SXP25.5 (Raioelimana na et al. 201 1) 

Tentative 

SXP565 (Bird et al. 2011) 

SXP293 (Schmidtk e et al. 2004) 

Tentative 

Possible new period 



oo 




Fig. 3. Color indicies of the candidate counterparts of the 
highly absorbed X-ray sources of our sample. For the sources 
reported with a filled circle, the brightest optical and infrared 
candidate counterparts are coincident, while the sources re- 
ported as an open circle correspond to those with the 'a' flag in 
Table [3] in this case, there is also a possible alternative coun- 
terpart, with a slightly brighter V magnitude, to the selected 
infrared counterpart. We also indicate the main-sequence stars 
and a bar which shows the reddening coiTesponding to Av=l 
mag. 



extractor), including background sky subtraction, of data 
allowed the creation of one-dimensional spectra, wavelength 
caUbrated using standard lines from a Cu-Ar arc. Observing 
conditions were not photometric, seeing was limited to 2" at 
best but varied throughout the evening. We show in Fig. |4] the 
resulting spectrum: its moderate signal-to-noise ratio does not 
allow us to make a spectral classification, but a significant Hq 
emission line, with equivalent width EW = -37.44 A, is clearly 
visible. 

Despite no periodicities being found in the OGLE light 
curve of this star, these findings indicate that source # 13 is a Be 
HMXRB. Since sources of this class typically have a transient 
nature, we searched for hard X-ray outbursts i n all the pub- 



lic data obtained with the INTEGRAL satellite (IWinkler et al 



2003|). that is all poi ntings obtained with the IBIS instrument 
(lUbertini et al.ll2003h in which the source position was within 
12 degrees from the center of the field of view. This led to a to- 
tal of 619 pointings, with different exposure times from about 



^i,f0404^^P4*li 



_J I L 

5000 6000 7000 

Waveleiigth (niigsLroms) 



'Vk 



Fig. 4. SAAO optical spectrum of source # 13; note the pres- 
ence of a strong H^ emission line. 

1800 sec to about 3600 sec, spanning from July 2003 to June 
2009. We used version 9.0 of the Off-line Scientific Analysis 
software to analyse the data. For each pointing we extracted an 
image in the 17^0 keV energy band, but in none of them was 
the source significantly detected. The typical flux upper limit 
of the individual images is on the order of ^20 mCrab. We also 
produced a deep image by summing all the individal pointings, 
reaching a total exposure of about 2 Msec. The source was not 
detected, with a 5 cr upper limit of about 1 mCrab in the 17-40 
keV range, which is well above the extrapolation of the flux 
seen by EPIC. 

6. Discussion 

We have detected 30 sources with hardness ratios indicating a 
column density Nh > 3 x lO^'^ cm~^. Five of these sources 
are HMXRBs, of which one was discovered in our search and 
spectroscopically identified as a new Be binary system. Most of 
the remaining 25 sources have very faint (V> 18) or no optical 
counterparts, with colors not compatible with reddened early 
type-stars (unless a significant IR excess is present). Therefore, 
they can hardly be considered as candidate absorbed HMXRBs. 
Spectral information for a few of the brightest sources of 
our sample was previously reported in other works. The XMM- 
Newton spectrum of the 25 s pulsar XMMU J0048 14. 1-73 1003 
was closely fit with a power law of photon index 1.33±0.27 
and Nh — 5 x 10^^ cm^^, confirming our conclusion that it 
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is a highly absorbed source. Our brightes t source (# 18) is the 
293 s pulsar RX J0058.2-723 1 , for which lHaberl et alJ (l2008ah 



reported an absorption N^ = (1.35 ± 0.32) x 10 cm~ , con- 
siderably smaller than the value we inferred from the hardness 
ratio. This might be due to the results of iHaberl et alj ( l2008al) 
being obtained from the June 2007 XMM-Newton observation, 
while our results for this source are based on the SMC Survey 
pointing carried out in October 2009, when the source was 
about a factor of ten brighter. HMXRBs can display large vari- 
ations in their intrinsic absorption, as shown for example in the 
case of the long period (1323 s) SMC pulsar RXJO 103. 6-7201 



(lEger & Haberi2008l) . 

We do not expect of course all the sources in our sample 
to be HMXRBs in the SMC. A large fraction of them will 
turn out to be background-absorbed AGNs. On the other, hand 
our sample should not contain foreground stars belonging to 
our Galaxy, which would show softer and less absorbed X-ray 
spectra. 

A rough estimate of the number of absorbed AGNs in our 
sample can be derived using published number flux relations 
for extragalactic sources. For example, based on the LogN- 
LogS rela t ions f or AGNs with N^ = 3 x 10^^ cm~^ by 
Gilli et al.l (120071) . we obtain a surface density of ^4 AGNs 
per square degree at a flux limit of 4 x 10^^"* erg cm^^ s^^. 
Our survey covered a sky area of ~^5 square degrees, although 
with a non-uniform sensitivity, because of telescope vignetting 
and different exposure times. Owing to the uncertainties in the 
LogN-LogS relations and the small number of objects, this pre- 
diction is consistent with the number of unidentified sources in 
our sample. 

7. Conclusions 

Motivated by the presence of several highly absorbed 
HMXRBs in our Galaxy, we have carried out a search for sim- 
ilar systems in the SMC exploiting the nearly complete cover- 
age of this galaxy obtained with XMM-Newton. Our selection 
criteria, corresponding to an absorption threshold of ^ 3 x 10^'^ 
cm^^, were met by four known HMXRBs and led to the dis- 
covery of a new highly absorbed Be binary in the SMC. We 
also selected other 25 additional sources, among which other 
highly absorbed HMXRBs might be present. Since we also ex- 
pect a good fraction of these sources to be AGNs, we can say 
that the SMC does not contain much more than ^10 persis- 
tent HMXRBs with intrinsic N/^ > 3 x 10^^ cm"^. We cannot 
exclude the presence of a much larger population of intrinsi- 
cally absorbed binaries of transient nature, which could not be 
detected in our survey if their quiescent luminosity is below a 
few 10'^* erg s^^. 
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